Entanglement generation in double-A system 



Ling Zhou, Yong Hong Ma, Xin Yu Zhao 
School of physics and optoelectronic technology, Dalian University of Technology, Dalian 116024, P-R-China 

In this paper, we study the generation of entanglement in a double-A system. Employing standard 
method of laser theory, we deduce the dynamic evolution equation of the two-mode field. We 
analyze the available entanglement criterion for double-A system and the condition of entanglement 
existence. Our results show that under proper parameters, the two-mode field can entangled and 
amlified. 
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I. INTRODUCTION 



Continuous variables entanglement (CVE), as entan- 
glement resource, has attracted lots of attention because 
CVE not only has advantages in quantum-information 
science [1] but also can be prepared unconditionally, 
whereas the preparation of discrete entanglement usu- 
ally relies on an event selection via coincidence measure- 
ments. Conventionally, continuous variables entangle- 
ment has been produced by nondegenerate parametric 
down- conversion (NPD) [2]. In order to improve the 
strength of the NPD, engineering the NPD Hamiltonian 
within cavity QED has also attracted much attention [3- 
5]. Besides parametric down- conversion [2, 6, 7], Xiong 
et. al. [8] had shown that two-photon correlated sponta- 
neous emission laser can work as a continuous variables 
entanglement producer and amplifier, which open a new 
attracting research domain. And then, a number of dif- 
ferent schemes have been proposed [9-14]. Different from 
the gain medium atoms in [8-14], Ref. [15] has studied a 
single-molecular-magnets system to produce CVE where 
physics process is similar to [11]. All of these works deal 
with the similar physics process where both of the two 
mode will be created (annihilated) a photon in one loop 
respectively ( similar to down-conversion system ). 

In this paper, we proposed a scheme to generate CVE 
where the one mode is created a photon and the other 
is annihilated, which is different from [8-16]. The system 
consists of atoms in double-A configuration interacting 
with two modes cavity fields. The atoms are driven into 
a coherent state of the upper two levels by two classical 
field. We obtain the master equation of the two mode 
fields. Through analysis of entanglement, we find that 
the criterion proposed in [18] can be used to judge en- 
tanglement. We show that in double-A system, entangle- 
ment exist on the condition that the two-mode quantum 
field is tuned away from the atomic transition, and the 
initial field is in a quantum state. Our study is help- 
ful to understand the entanglement charateristic within 
a system where quantum field is in "y configuration. 



II. THE MODEL AND THEORY 
CALCULATION 

We consider a system of atoms in double-A configu- 
ration shown in Fig.l. Two cavity fields interact with 
atomic transition \a) |c) and |6} <-> |c) with detun- 
ing Aq and Af), respectively. The two classical pumping 
fields with Rabi frequency 1^2 and f2i drive the atomic 
level between \a) ^ \d) and \b) ^ \d) with detuning 
Al and A2 respectively. Our double-A system can be 
sodium atoms in a vapor cell [28] where the lower states 
are the two hyperfine levels \F = 1) and F = 2) of 3^S'i/2, 
and the upper state are \F = 1) and F = 2) of 3'^Pi/2- 
The double-A system also can be atomic Pb vapor [19]. 
The phase-dependent electromagnetically induced trans- 
parency [28] and efficient nonlinear frequency conversion 
[19] have been investigated experimentally in double-A 
system. Ref. [20] studied dark-state polaritons in double- 
A system. Here, we are interested in producing two-mode 
entangled laser via the double-A system. In interaction 
picture, the Hamiltonian of the system can be written as 



Hq = iyia[ai + 1/20, 



2^2 



(1) 



-|-a;o|a}(a| + i^b\b){b\ + ujc\c){c\ + uJd\d){d\ 



Hi = giai\b) {c\ + g2a2\a) {c\ 

+ni\b){d\e~''^'' + n2\a)(d\e- 



(2) 



H.c. 



We hope that the Hamiltonian do not contain time t so 
as to simplify the density matrix deduction of the field. 
In order to do that, we assume that the classical fields 
detuning Ai = Aa — A and A2 ~ A^ —A. Now we 
goes into a frame by performing a unitary transformation 
U = exp{i[Ho + Aa|a)(a| + Ab\b){b\ + A\d){d\]t}. In the 
new frame, the Hamiltonian is as 



Hi = -Aa|a)(a| - A6|6)(6| - A|d)(d| 
+ [5iai|6)(c| + g2a2\a){c\ 
+ni\b){d\e-'"^^' + n2\a){d\e-'"^^ 



(3) 



H.c.]. 



In order to see the entanglement of the two-mode field, 
we need to obtain the equation of motion of the two- 
mode field. Using the standard procedure in laser theory 
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and 



c) 



FIG. 1: The level configuration of atoms. Two cavity modes 
interact with atomic transition \b) ^ |c) and \a) ^ |c) with 
detuning Ai, and Aa respectively while the two classical fields 
drive the atomic level between \b) ^ \d} and \a) ^ \d) 
with detuning A2 and Ai. For simplicity, we assume the 
spontaneous-emission rate of four level are the same. 



developed by Scully and Zubairy[16, 21, 22], we obtain 
the following master equation governing the dynamics of 
the two-mode cavity fields as 

p = -Ki{a\aip ~ aipa\) - K2{ala2p - a2pal) (4) 
—ai{paia\ - a\pai) — a2{pa2a\ — alpa2) 
-ai2{pa2a\ - a\pa2) - a2i{paial - alpai) + h.c. 

We can see that the master equation has the term 
pa2a\—a\pa2 which means that the one mode is created a 
photon and the other mode is annihilated a photon. The 
detail deduction of the equation is given in appendix A. 
In Eq. (4) , we have include the loss of the two- mode cavity 
with loss rate ki and K2- The coefficients are 



ai 

a2 
ai2 

0:21 



9i 

— [AliLbb -I- A2lLab + A^iLdb], 

^[Ai2Lba + A22Laa + 
glg2 

D 
.91.92 

D 



(5) 



[AiiLba + A2lLaa + A^iLda], 
[Ai2Lbb + A22Lab + A^2Ldb\, 



where Z? = (7 — i/S.a){l ~ i^b)[l ~ + 
^l{-f - iAa) + nlij - iAb), and 
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Laa 
Lab 
Ldb 
Lda 



bb 



— ii^2 

1 

7 + ^(A2 - Ai) 



7 



27 
7 - »A2 

27 
7 - *Ai 



27 



y2 



-yi 

iQ2 



-2/2, 

iQ2 
27 



(6) 



2/2 



-2/3, 



27 



2/1, 
2/1, 



All 
A12 
A22 
A32 



A21 ^ -nin2, A31 = -ii7i(7 - iAa), 

(7-^A,)[7-^(A+^ + ^], 
-^^12(7 ~ *Ab) 



(7) 



with 



2irins{aifl2 - bfli) 2irins{a2^i - 6^^2L^ 
2/2 = — r9 '2/3 = — u (8) 



aia2 



-52 



aia2 



62 



2/1 = 
in which 



r!2(Ai-2A2) rii(2Ai-A2) 
2/2 + 2/3, 



s 


= 72 + 172 + ^2 + 


(A2-Aif, 


ai 


= Mis-fll{2A2 


-Al)^ 


a2 


= Af2S-r22(2Ai 


-A2)2, 


b 


= r^if^zps - (Ai 


-2A2)(2Ai-A2)] 


Ml 


= 72 -1- 4^2 -1- f72 


+ A2, 


M2 




+ A2. 



(9) 



Although our four- level atom is similar to [11, 15], the 
physical process of the two-mode quantum fields is dif- 
ferent because the two quantum fields work in different 
atomic level. In [8-15], both of the two mode will be cre- 
ated or annihilated a photon in one loop. So the master 
equation is of the form pajoj — a\pa2 (^02^1 — o,ipO'2)- 
In our system, the two quantum fields are in a " V" form 
levels if we do not see the two classical pumping fields. 
The simplified " V" form levels is similar to " Hanle effect" 
laser [? ] where the master equation is with the term 
pa2a\ — a\pa2. In our system, the two classical fields 
make the atoms with the coherence of the two up-level 
\a) and \b) [see (AlO)]. When the spontaneous emissions 
from \a) and \b) to |c) take place, entangled photons will 
be produced. 



III. ENTANGLEMENT CRITERION CHOICE 
AND THE DISCUSSION OF THE 
ENTANGLEMENT CONDITION 



How to determine the entanglement is a key problem. 
In Ref.[8-15], employing the criterion (Au)2 + (Aw)2 < 
2 [23], a inequality of the sum of the quantum fiuctua- 
tions of two operators u and v for some entangled state, 
they find the entanglement between the two mode fields. 
However, the criterion inequality of the sum of the quan- 
tum fluctuations can not be applied to measure coherent 
state[24]. Although the entanglement criterion on mea- 
sure continuous variable have been developed[22-25], we 
still can not find a criterion to judge all kind of contin- 
uous variable entanglement. In order to make clear the 
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kind of entanglement existing in our model, wc now dis- 
cuss the analytic solution in our system so as to choose 
a appropriate entanglement criterion as well as to know 
the condition of entanglement. 

Now we analyze the entanglement condition . If = 
(ji2, ^1 = ^2, and A = Aa = Afe » ^i,fl2,^, through 
Eq.(6) to (10), one can obtain the relation ai = a2 = 
= ct2i = ict {a is a real number). Usually, the loss 
of the cavity do not change the entanglement structure 
of the state. It just destroy or sometimes enhance the 
entanglement a little. So, in our choice entanglement 
criterion, we omit the loss of the cavity. Therefore, the 
master equation of our system Eq.(4) can be simplified 
as 

p = ia[aia\ + a2a\ + a2a[ + aia\^ p\. (10) 

The effective Hamiltonian H^q — —a{aia\+a2a\+a2a\ + 
aiaj). Due to [aia\+a2a\,a2a\+aia2\ = 0, in interaction 
picture Hgqi = —a{a2a\+aia2)- One can easy check that 
the system state, evolved by H^gi = — a(a2a| + aiOj), 
never meet with the criterion (A?i)^ + (Af)^ < + -p- 
for the initial field number [ni,ri2). We recognize the field 
Hamiltonian is the generator of the SU (2) coherent state 
[27]. The evolution of the state |*(0)) is 

\-qi{t)) = e-*^"^*|*(0)) = e^+-^+e^«^'^^«e^-^-|*(0)), 

where = a\a2, = aia\. These operators sat- 
isfy the SU{2) commutation relations, i.e., [iiT- , = 
-2Ko, [Ko,K+] = K+, [Ko,K_] = -K_, with Kq = 
^{a\ai — ala2); and in which 

xq = {cosh iat}~ ^ , 
x+ = X- = tanh iat. 

If the initial field state is two-mode Fock state |0, A''), the 
evolution of the state is 

|4'(f)) = {cosatf/'^ Y.[ n I (*tanQt)"|n, TV - n). 

71=0 ^ ^ 

(11) 

From the entanglement definition of pure state, we know 
that the state \^{t)) is a entangled one. 

Unfortunately, the Hamiltonian Hgqi can not entangle 
initial coherent state, because the evolution of the system 
as 

l^-(t)) =e-+^+e^'''"-''e--^-|/?i,/32> = |^i,^2>, (12) 
with 

Pi = cosat -|- i/32 sinaf, 
02 = 1^2 cos at + ipi sin at. 

So, it is not entangled. 

The two-mode SU{2) cat state is sub-Poissonian dis- 
tribution. We recall the criterion, proposed by Hillery 



and Zubairy[18] can be used for non- Gaussionian state. 
The criterion say that if 

{N1N2) < |(ai4)|2 (13) 

the two-mode field is entangled. If the field initially is 
in number state |ni,n2), using the differential equation 
Eq.(Bl-B13)(let «; = and ai = a2 = Q!i2 = ^21 = ia), 
we finally obtain 

{N1N2) — 1(0102)!^ = nin2 — ^{ni+n2 + 2nin2) sin^ 2at. 

(14) 

The maximum value of sin 2at is 1 ; therefore if 

2nin2 <ni+ n2, (15) 

the two mode field will be entangled. Because ni and 
n2 are integer, in order to meet with 2n\n2 < rii -|- n2, 
the number m and n2 should be not equal. If either 
ni or n2 is zero (the state is standard SU{2) coherent 
state), we can see that (NaNt) — \{ab^)\'^ is always less 
than zero; thus we say the state is entangled. Therefore, 
the criterion Eq.(12) can be used for judge entanglement 
within our system. 

However, for resonant case (A5 = Aa = A = 0), 
if 7h = 7a ' 51 = 52 and fli = Q2, the coefficients 
ai = a2 ~ a\2 — 0121 — f3 (real number). For the 
initial state \ni,n2), after complicated calculation em- 
ploying Eqs. B1-B13 for ai = 02 = ai2 = a2i = P, we 
have 

{NiN2)-\{aial)f = 2l^i3e»^^ + 2e^^' ~ 

8 

+ ^{l + e^^*-2e^^*)h0 

If rii or 712 is zero, {N1N2) — \ {aial)\'^ equal to zero at 
initial time. Except that the {N1N2) — \{aial)\'^ is larger 
than zero. It is obvious that we can not obtain entan- 
glement in resonant case. This conclusion is consistent 
with the work in Ref.[12], where author show that for 
two-level qiiantum beat laser, entanglement can be cre- 
ated only when the strong driving field should be tuned 
away from the atomic transition. 

IV. THE ENTANGLEMENT OF THE CAVITY 
FIELD 

In above section, we discuss a special case so as to 
choose entanglement criterion and make clear the condi- 
tion of entanglement. Although above analysis is for pure 
state (approximation of master equation Eq.(4)), But the 
criterion {N1N2) < |(aia2)P should be available in judg- 
ing entanglement for general case. Now, considering the 
loss of the cavity and the decay of the atomic levels, we 
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FIG. 2: (a): The time evolution of the entanglement. (b): The 
time evolution of the two-mode fields where red line is for A'^i 
and blue line is for N2. Initially, the atom is in number state 
|10, 0) The parameters are gi = ff2 = 1; Aa = A;, — 50, 
A = 4, n„ = 20, 7 = 1, Ki = K2 = 0.010. = S^i = 5. 

numerical solve the differential Eqs. (Bl) to (B13) and 
plot the entanglement criterion {N1N2) — Haia^)]"^ and 
the Ni{{a\ai)),N2i{ala2)). 

In Fig. 2, we plot the case that the initial field state is 
a number state |10, 0) where Af, = Aa » 7{, = 7^ which 
means that the classical field resonantly drive the atom 
(Ai = A2 = 0) and the quantum field interact with the 
atoms with equal detunings. We see that due to the loss 
of the cavity, the entanglement gradually disappear and 
photon number of the two-mode field also decrease under 
large detuning case. Of course, if the cavity is ideal, one 
will observe the entanglement oscillation. 




FIG. 3: The time evolution of entanglement. Initially, the 
atom is in number state |1, 0) The parameters are gi = g2 ~ 1; 
A„ = 50, At = 20, A = 10, nn = 20, 7 = 1, ki = k2 = 
0.01. Q2 = ill = 4,5,6 for dotted, dashed and solid line, 
respectively. 

However, with the same detuning Ah = Aq (when 
<?i — 92), we can not have amplified entangled laser 
shown in Fig. 2. The quantum fields are in "V" form. 
If Aft = Aa, the photon number in two mode only oscil- 
late because of the symmetry. In our numerical simula- 
tion, we find that in order to have amplified entangled 



laser, Aq and A;, should be different. For initial field 
state in number state |1,0), we plot entanglement and 
average photon numbers in Fig. 3 and 4 for several val- 
ues of ill ($12). One can see clearly that entanglement can 
be obtained without preparation atomic coherence before 
(here, atoms are injected in state \d)). But the photon 
number in two mode has large difference. By adjusting 
the values of fli(fl2), we can adjust the time region of en- 
tanglement. Because we inject the atom in atomic state 
\d), it will need time to evolve into a coherence among 
the atomic level |a), \b) and \d). So, we have no entan- 
glement during a initial short time . With large value 
of f2i(f22), the atoms will acquire their coherence quickly 
so that the entanglement appear quickly. However, with 
large value of ili{n2), the photon number also will be 
amplified quickly shown in Fig. 4. As in our analytic 
calculation, we have known that the photon number in 
two mode differ (Eq.(15)). Here, in order to amplify the 
photon number, the photon number not only should have 
difference but also can not put up with very large pho- 
ton number. With the increasing of photon number, the 
entanglement disappear. But the disentanglement is not 
resulted from loss of the cavity because we find even for 
K = 0, entanglement also disappear. We conclude that 
the disentanglement result from the increase of photon 
number rather than from the loss of the cavity. As it is 
pointed out in Rcf. [21],in the high-gain limit the condi- 
tion in Eq.(9) is no longer able to detect whether there 
is entanglement in the state. 




FIG. 4: The time evolution of average photon number. 

Now, we show another function of the classical 
fields, i.e., the ability to overcome the loss of the cavity 
which is shown in Fig. 4. Let us compare dotted line and 
solid line. The two lines correspond to the loss rate of 
the cavity ki = K2 = 0.01 and 0.1, respectively; and all 
the other parameters are the same. Due to the increasing 
of ki(k2), the values of {N1N2) — |(aia2)P move up. If 
ki(k2) keep increasing, we will loss entanglement. How- 
ever, with the help of classical fields, we still can obtain 
entanglement even through ki(k2) is large, which can be 
observed by comparing dashed line and solid one. Al- 
though the loss rate ki = K2 = 0.1, through increasing 
ili{n2) to 6, we still can have entanglement. Of course. 
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because of the increasing of ^li{^2), the time region move 
left, which we have analyze it in Fig. 3. 



APPENDIX A: CALCULATION DETAILS OF 
DENSITY MATRIX OF TWO-MODE FIELDS 




FIG. 5: The time evolution of entanglement. Initially, the 
atom is in number state jlO, 0) The parameters are <7i = (?2 = 
1; Aa = 50, At = 20, A = 10, n„ = 30, 7 = 1. Dotted line: 
Ki = K2 = 0.010. Q2 ~ ^1 ~ 4;solid line: ki — k,2 = 0.1, 
i}2 = i^i — 4; dashed line:Ki — K2 = 0.1, Q2 ~ ^1 = 6. 



CONCLUSION 



In conclusion, we have studied the generation of entan- 
glement in a double-A system. We derive the theory of 
this system and analyze the available entanglement crite- 
rion for double-A system. When the atoms are injected in 
the ground state \d), the entangled laser can be achieved 
under the condition of suitable parameters. Due to the 
classical pumping field introduction, we do not need to 
prepare atomic coherence, and the intensity of the quan- 
tum fields will be amplified. The classical pumping can 
overcome the loss of the cavity. Our results show that the 
time for which the two modes remain entangled depends 
upon the strength of the Rabi frequency of the classical 
driving field. 

Our results is helpful in understanding the entangle- 
ment characteristic when the master equation contain 
the term pa2a\ — a\pa2 such as quantum beats laser and 
Hanle effect laser system. Our studies is limited to the 
initial state [1,0). One can research other initial field 
state. Our initial field should be easy to obtain. Let ex- 
cited two-level atom with transition frequency vl{ or v2) 
passing through the vacuum two-mode cavity, when we 
detect the output atom in ground state, we will have the 
field state |1,0). 

Acknowledgments: Authors thank Professor M. S. 
Zubairy and M. Ikram for their critical reading. 
The project was supported by NSFC under Grant 
No.10774020, and also supported by SRF for ROCS, 
SEM. 



The classical fields will be treated to all orders in the 
Rabi frequency. The transitions \a) - |c) and \h) - \c) 
are treated fully quantum mechanically but only up to 
second order in the corresponding coupling constants. By 
partially tracing the global state of Schrodinger equation 
over the atomic variables, we have the formal reduced 
fields 

Pf = -^{[Hcb,Pbc] + [Hca, Pac] + Wac-, Pea] + Wbc-, Pcb])- 

(Al) 

Eq.(Al) reveals that we need to get the density matrix 
elements , etc.. Inserting the Hamiltonian Eq.(l) 
into A.l, from the schrodinger equation, we have 

Pbc = -{l~i^b)Pbc-i^iPdc (A2) 
+ {-'igiaiPcc + igiPbh^i + ig2Pbaa2), 

Pac = -~il-i^a)Pac~i^2Pdc 

+ {i9lPabO'l + i92Paaa2 " ig2a2Pcc), 

Pdc = -[l-i^]Pdc-i^lPbc~i^2Pac 
+ {mPdbai+ig2Pda'^2)- 

In the last equations Eq.(A2), we have consider the 
spontaneous-emission of the atomic level. We rewrite it 
in a matrix form as 



P 



-Mp + A 



where 




M 



— iAf, iili 
7 - iAa i^2 



WlPbb^l + tg2Pbaa2 
^ = I mPabO-l + ig2Paa0'2 
"i-giPdbO-l + i92PdaO-2 



(A3) 



(A4) 



(A5) 



(A6) 



When we write matrix A, we let p^^ = and will explain 
the reason later. A solution of Eq.(A3) which is a linear 
in the coupling constant gi(^2) can be obtained [16, 21, 22]. 
Here we only care for the matrix elements and p^^,, 
so we just write the solution of the two terms as 

Pbc = -^[{AiiPbb + ^21 Pab + ^31 Pdb)gia'i 

H^llPba + Mipla + ^3lPda)52a2], (A7) 



Pac ^ -^[{Al2plb+ A22plb + A^2P^db)gia.l 

+ {Al2pla + A22pla + ^32pl)ff2a2] (A8) 
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with 



All 

Az2 



(7-iA„)(7-iA) + 02, 

A21 = -flifl2,Asi = -irJi(7 • 
(7-iAb)(7-iA) 
-ifl2ij - i^b), 



(A9) 



where = (7 - «Aa)(7 - zAfa)(7 - iA) + (7 - iA^) + 
^2(7 ~ ^'Afc) ill Eqs.(A7)(A8). As a approximation, the 
density matrix elements in right side of Eqs.(A7)(A8) 
such as p"^, p^^, etc. will be determined by steady state 
of classical fields. In other words, the density matrix ele- 
ments pIJj,, p"^, etc. of classical fields, as a zero order ap- 
proximation, are substituted into right side of Eqs.(A7) 
and (A8) . And then, we can obtain a first order approx- 
imation of density matrix elements , in terms of 
couphngs gi (32) ■ 

Now, we just consider classical fields to determine the 
zero order approximation of the density matrix elements 
Pbb^ Pba^ '^^^ differential equations of density matrix 
elements only with classical fields and atomic decay are 



-7pg,-ini(pS(,-pg<,), 

-lPaa-i^2{Pda- Pad), 

-[7 - z(A2 - A,)]pl + i^2Pld - ^f^iPl, 

-[7 + *Ai)]p1 - in2{pla - pIi) - ^^iPla, 

-[7 + iA2)]p°, - zf^2P°i, - i^hiplb - PmI 



Pbb = -IPbb-i^Mb-Pbd), (AlO) 

Paa 

■pt ^ 
p1 

Pdd = -lPdb-i^l{Pbd-Pdb)-i^2{p°ad-Pda) + nnP 

with Ai = Aq — A,A2 = Aft — A. For p^^, we have 
p"^ = —^Pcc- The steady state solution p^^ = (It is 
the reason why we let p^^ = in Eq.(A6)). Substituting 
the steady state solution of (AlO) into (A7) and (AS), 
we obtain p^^ , p^^. And then, inserting p^^ and p^^ back 
into (Al), one can have the master equation Eq.(4) with 
coefficients Eq.(5) to (9). 



APPENDIX B: CALCULATION DETAILS OF 
DENSITY MATRIX OF TWO-MODE FIELDS 

In order to numerical calculate the entanglement cri- 
terion {N1N2) - |(ai4)p and the Ni{N2), we need to 
deduce a series differential equations from master equa- 
tions (4) which are listed below. 



' = {ai+al^2Ki){a\ai) (Bl) 

+Q;i2(«2"i) + ai2(a2a|) -|- ai -|- aj, 

= {a,+a*2-K,-K2){a,al) (B2) 

^ ^ +a2i(aiai) -I-Q!i2(4a2) -|-Q!i2 -|-a2i, 

- = \fiL\-\- a2-\-2(y.2 — i^\ — i>>^2)\a2^2<Ho-\) 



dt 



dt 



d(aia^aia^) 
di 



+{a2 + ^2 + 2K;2)(a2^ai) -|- 2a2i(a2a2ai'ai) 
-|-2a2i(a2'^2) + Q!2i(a2^of) -|- ai2(a2a2a2a2^) 
+a*2i{a+ai)+a*2^, (B3) 



2{al + a2-K2- Ki)(a^a+2) (B4) 

-|-2a2i(alaia^"a2) -I- 2(a2i + Q;i2)(a^a2) 
-|-2a^2(<^2a2'^2af), 

= 2(q!i +al- 2Ki){aiataiat) (B5) 
+{ai +al+ 6Ki)(afai) 
+2ai2{a\aia^ a2) + 20*2(010^0102) 
+ai2{a2ai) + 0:^2(0102) 
-|-4ki + 0:1+0;^, 

(ai + al+a2+a2- 2ki - 2K2){afaiala2) 

+{a-i + a*)(o^a2) -|- (a2 + a2){alai) 
+(X2i{aiaia\a2) + a2i(oiaf oioj) 

+Q!l2(a202'^20f)Q!l2(a2'^202'0i). (B6) 



Substituting the subscript 1 (2) with 2 (1) and then mak- 
ing their Hcrinitian conjugate through (Bl) to (B5), we 
can obtain the other seven differential equations. The 
totall thirteen differenttial equations will be a closed set. 
We can numerical solve it. 



^(0^010202) 



dt 
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